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ABSTRACT 

We discuss the properties of spiral arms in a N-body simulation of a barred galaxy and present 
evidence that these are manifold-driven. The strongest evidence comes from following the 
trajectories of individual particles. Indeed, these move along the arms while spreading out 
a little. In the neighbourhood of the Lagrangian points they follow a variety of paths, as 
expected by manifold-driven trajectories. Further evidence comes from the properties of the 
arms themselves, such as their shape and growth pattern. The shape of the manifold arms 
changes considerably with time, as expected from the changes in the bar strength and pattern 
speed. In particular, the radial extent of the arms increases with time, thus bringing about a 
considerable increase of the disc size, by as much as 50% in about a Gyr. 

Key words: galaxies: spiral - galaxies: structure - galaxies: kinematics and dynamics - 
galaxies: evolution. 



1 INTRODUCTION 



In a series of papers i Romero-Gomez et al. |2006[ Paper 



Paper IV; Athanassou 



I; Romero-Gomez et al. 2007 Paper II; [Athanassoula, Romero 



Gomez & Masdemont 2009a Paper III; Athanassoula et al. 2009b 



ro- 

3b 



a et al.|2010[ Paper V), we proposed a the- 



ory to explain the formation and properties of spirals and inner and 
outer rings in barred galaxies. According to it, the backbone of 
these structures are a bunch of orbits guided and confined by the 
invariant manifolds associated with the periodic orbits around the 
saddle points of the potential in the frame of reference co-rotating 
with the bar. We call our theory manifold theory, or manifold flux- 
tube theory. 

Some of the introductory dynamics necessary to follow our 
work is summarised in (Binney & Tremaine 2008, Sect. 3.3.2). 
Manifolds and the orbits they guide are described and explained 
in Paper I, while a relatively lengthy summary, avoiding equations, 
can be found in Sect. 2 of Paper III. Here we analyse a N-body 
simulation of an evolving barred galaxy and present evidence that 
its spiral arms are manifold-driven. In Sect. [2] we give a brief the- 
oretical reminder and describe relevant manifold shapes in simple 
analytical potentials. In Sect.[3]and|4]we present the simulation and 
our results and make comparisons with our theoretical predictions. 
Further discussion and conclusions are given in Sect. [5] 



2 THEORETICAL REMINDERS AND EXTENSIONS 

Let us model a barred galaxy potential in the simplest possible 
way, namely by an axisymmetric part (including the disc and the 
halo) and a bar, which we assume rotates with a constant angular 
velocity. The dynamics of this system is best studied in a frame 
co-rotating with the bar, where there are five equilibrium points at 



which the derivative of the potential in the rotating frame is zero. 
These are called Lagrangian points (Li, i = 1,5). Two of them (Li 
and L2) are located on the direction of the bar major axis, outside 
the bar but near its ends. By convention, we place the bar along the 
x axis, measure the angles in a counterclockwise sense and call Li 
(L2) the Lagrangian point on the right (left) of the centre. Li and 
L2 are saddle-point unstable, so that the families of periodic orbits 
around each one of them (called Lyapunov orbits, Lya punov|1949| > 
are, at least in their neighbourhood, also unstable. They thus can 
not trap regular quasi-periodic orbits around them since any orbit 
with initial conditions in their immediate vicinity (in phase space) 
is chaotic and will have to escape their neighbourhood. Not all de- 
parture directions are, however, possible. The direction in which the 
orbit can escape is set by what are called the invariant manifolds. 

Manifolds can be thought of as tubes that guide the motion 
of particles whose energy is equal to theirs. There are four man- 
ifold branches emanating from a given Lyapunov orbit (Fig. 2 of 
Paper III), two inside corotation (inner branches) and two outside 
(outer branches). Along two of these branches (one inner and one 
outer) the mean motion is towards the region of the Lagrangian 
point, while along the other two it is away from it. In Papers I - V 
we proposed that these manifolds and the orbits they guide are the 
building blocks of the spirals and rings in barred galaxies. 

Fig.[T|illustrates a number of useful manifold shapes and prop- 
erties. In the left panel, we use a potential including a spiral com- 
ponent. The latter exerts an additional gravitational forcing which 
is not symmetric with respect to the bar major axis, thus shifting 
the Lagrangian points away from that axis. The motion along the 
inner manifold branch (light green), is from the upper bar region 
to the left, i.e. towards the L2, as shown by the black arrow. These 
manifolds then circle for a few times around the L2 and then fol- 
low the outer manifold branch (red), in the sense shown by the ar- 
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Figure 1. Manifolds related to the L2 Lagrangian point in two simple analytical potentials. We plot in light green the inner manifold branch, in red the outer 
one and in blue a specific manifold (see text). In black we plot the zero velocity curve, and we mark by arrows the direction of the motion along the manifolds 
and with asterisks the positions of the Li, L2 and the galactic centre. The thin black solid line indicates the direction of the bar major axis and the dashed one 
the direction joining the Li to the L2. 



row, thereby outlining the spiral. After tracing somewhat more than 
180° in azimuth, they form a loop each, as if they were bouncing 
off a circle with radius roughly equal to the outer distance of the 
zero velocity curve from the galactic centre. All displayed mani- 
folds have the same Jacobi constant, near that of the Li and L2 . 

In the middle and right panels, we use a simple bar potential 
with no spiral component, so that the Li and L2 are now on the 
direction of the bar major axis. The manifolds are very similar to 
those of the previous example, with one considerable difference, 
i.e. we chose now to display manifolds whose Jacobi constant is 
much larger than that of the Li and L2. Therefore, both the Lya- 
punov orbits and the outline of the corresponding manifolds are 
much larger than in the previous example (see also Paper I) and 
part of the outer (red) manifolds, when reaching the vicinity of Li 
will turn inwards, becoming inner manifolds and moving leftwards 
along the upper part of the bar towards the L2 . We plot, in blue, one 
of these manifolds as an illustration. Note that amongst all the red 
manifolds, the one that turns inwards is the one which is leftmost 
in the last part before reaching Li, This motion is characteristic of 
manifolds with relatively large Jacobi constant and is explained in 
detail in K oon et aTj pOOO) and in Paper II (cf. also Paper V). The 
fraction of orbits that follows such trajectories depends on the po- 
tential and the distribution function of the orbits. In the right panel 
we follow the evolution of the blue manifold further to later times 
and show that it retraces the upper outer manifold branch. Another 
alternative (not shown here) is that, in the neighbourhood of the L2, 
instead of going upwards and following the arm, the trajectory goes 
downwards and follows an inner branch below the bar major axis 
(as e.g. in Fig. 1 of Paper 5, or the black path in the left panel of 
Fig. 3 in that paper, albeit for a different potential). 




Figure 2. Model rotation curve (thick solid line). The halo and disc contri- 
butions are given by a dashed and a thin solid curve, respectively. 



where r is the radius, Mh is the halo mass, 7 and r c are the halo 
core and cut-off radii, respectively, and the constant a is given by 
a = [1 — \pK ex P(q 2 ) (1 — er f{l))]~ > where q = 7/r c 
(Hernquist 1993 ). We use here 7 = 15 kpc, r c = 50 kpc and 
M h = 25 x 10 1 " M . The initial density distribution of the disc is 



Pd(R, z) = — ^ exp(~R/h) sech 2 ( — ) 

VKtl £ Zq Zq 



(1) 



where R is the cylindrical radius, M4 is the disc mass, h is the disc 
radial scale length and zo is the disc vertical scale thickness. We 
use here h = 3 kpc, z = 0.6 kpc and M d = 5 x 10 10 M . 
The corresponding rotation curve is shown in Fig. [2] For the radial 
velocity dispersion of the disc particles, <jr(R), we take ur(R) — 
100 • exp(-R/3h) kms" 1 . 

The initial conditions were made using the iterative method 
(Rodionov, Athanassoula & Sotnikova 2009), and the simulation 
was run using the GADGET2 code (Springel, Yosh ida, & White] 
[200T]|Springel|2005> . We adopted softening lengths of 100 (200) 
pc for the disc (halo) and an opening angle of 0.5. 



3 SIMULATION 

The simulation we will discuss here has two live components, a 
halo and a disc, represented by a million and 200 000 particles, 
respectively. The former has a volume density 

( ) - Mh 2L ex P(- r2 / r c) 

Ph{r > ~ 27rV 2 r c r 2 + 7 2 ' 



4 RESULTS 

To follow the morphological evolution in this simulation, we saved 
the snapshots every 0.005 Gyrs. For each one of them, the position 
angle of the bar was calculated and the snapshot rotated so as to dis- 
play the bar horizontally. In this manner, we visualise the evolution 
in a frame of reference co-rotating with the bar. An animation, us- 
ing these frames and produced with the GLNEM02 software, can 
be viewed in |http://195.221.212.246:4780/dynam/movie/MFolds1 
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Figure 3. Evolution of the disc component. The time in Gyrs is given in the bottom left comer of each panel. 



sgs058_noorbits.avi and sgs058_polar.avi for face-on Cartesian and 
polar views, respectively, while selected frames are shown in Fig. [3] 

For our initial conditions, and within a distance of the order 
of 20 kpc from the centre, the dynamics are dominated by the 
disc rather than the halo mass distribution. Because of this, the bar 
grows early-on in the simulation and fast in time. It is initially very 
fat, but by t = 0.5, it starts becoming longer and narrower, while 
two spiral arms start forming, one from each end of the bar. These 
spirals are trailing and their angular extent increases with time as 
their tip approaches the opposite side of the bar. Until about 0.9 
or 1 .0 Gyrs they form a grand-design two armed structure, staying 
attached to the end of the bar and rotating with the same pattern 
speed as the bar. 

Shortly after t = 0.9 Gyrs, the shape of the bar undergoes 
strong changes - as material initially in the arms is accreted to its 
outer parts - while the arm-interarm density contrast drops, so that 
the spiral is not as clearly discernible. This is followed by a second 
spiral episode, qualitatively very similar to the first one and start- 
ing between 1.0 and 1.1 Gyrs. I.e. a new two-armed grand-design 
spiral develops, again starting from the tip of the bar and extending 
towards the opposite side of the bar. This second spiral episode lasts 
till about £=1.5 Gyrs and its shape and amplitude are quantitatively 
considerably different from those in the first episode. 

The above description is in good agreement with our mani- 
fold theory results. Indeed, this theory leads naturally to two-armed 
trailing spiral arms (Paper IV), which start growing from the tip of 
the bar first outwards and then towards the opposite end of the bar 
(Fig. 6, Paper I). Material for these arms should come from the 
outer parts of the bar (Paper V). One of the extremities of each arm 
should be linked to L2 or L2, and the arms should rotate with the 
same pattern speed as the bar. All these developments and morpho- 
logical properties are indeed seen during the evolution in our sim- 
ulation, arguing strongly for a manifold origin of the spiral arms. 



The main argument, however, in favour of the manifold origin 
of the spirals comes from the orbits of the individual particles. In 
density wave theory, the arms are loci of density maxima. Particles 
should thus traverse the arms, but stay longer in the arm than in 
the interarm region {Lin & S hu 1964). This is totally different from 
our manifold theory, where spiral arms should be a bundle of orbits 
guided by the manifolds, so that particles should move along the 
arms rather than across them. 

It should thus be possible to find out which of the two theo- 
ries is the main driver of the spiral structure in the simulation sim- 
ply by following a number of particles. For this we use our series 
of snapshots in which the bar orientation is kept horizontal. At t 
= 0.8 Gyrs, when the bar and spiral are well developed, we se- 
lected 60 particles located in the part of the spiral arm which is 
near the tip of the bar, roughly where we estimated by eye that 
the L2 Lagrangian point would lie (cf. left panel of Fig. [TJ. Sim- 
ilar results can be found by selecting particles at other times or 
locations, provided they are clearly in an arm at selection time. 
We then followed the trajectories of these particles from t = 0.5 
to 1.55 Gyrs and produced a sequence of 211 frames in each of 
which we superposed on the snapshot of all disc particles a filled 
white circle marking the current location of each chosen parti- 
cle and a white line for its trajectory over the previous 0.175 
GyrQ We repeated this task for all snapshots and thus produced 
an animation jhttp://195.221.212.246:4780/dynam/mo vie/MFolds 
sgs058_orbits.avi and sgs058_polar_orbits.avi for Cartesian and po- 
lar views, respectively). In Fig. [4] we show nine such frames which 
display the salient features of the particle trajectories. 

Roughly from t = 0.5 to 0.65 Gyrs the chosen particles travel 

1 This time range should be kept in mind when comparing the spiral shape 
and the early parts of orbit, because the spiral could have evolved during 
these 0.175 Gyrs. 
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Figure 4. Nine snapshots of the disc component, on which we overlay the locations (white filled circles) and trajectories (white solid lines) of the 60 particles 
that we follow (see text). The time in Gyrs is given in the bottom left corner of each panel. 



along the outer part of the bar, in the direction from the L\ to the 
L2, then they circle a couple of times around what should be the 
Z/2 (frame at t = 0.82) and escape its vicinity following the arm 
(times 0.90 to 1.18). Thus the particle orbits stay within the arms, 
as expected for the manifold theory, and do not cross them as would 
have been expected by the density wave theory. After roughly t = 
1.1 Gyrs, the forerunners of the group of particles we follow reach 
the vicinity of the opposite side of the bar from which they em- 
anated (i.e. the right side, near Li). From that point onward the 
particles divide themselves into two groups. Some continue down- 
wards, make a loop, as if they were bouncing off an invisible bar- 
rier, and then trace the lower spiral arm from the L\ to the L2. 
Others retrace the upper outline of the bar leftwards in the direc- 
tion from Li to 1/2. Once they have reached the vicinity of L2 an- 
other branching occurs. Some particles go to the upper spiral arm 
which they follow from the L2 to the L\, while others stay in the 
bar (frames at t = 1.45 and 1.53). All the trajectories and branch- 



ings described above are very similar to those seen in Fig.[T|and in 
Papers I - V. 

Furthermore, note that the orbits spread out as they trace the 
arm moving from the L2 to the L\ between t = 0.65 and 1.1 Gyrs, 
leading to a gradual widening of the spiral arm. This can be already 
seen in Fig. [3] This widening also is expected from the manifold 
origin of these arms (Fig. [T] here and Papers I and V). How much 
this widening is depends on the potential, so that we can not com- 
pare Figs.[T|and|4]quantitatively. 



5 DISCUSSION AND CONCLUSIONS 

In this letter we followed the formation of spiral structure in a 
N-body simulation of a barred disc galaxy. We witness the forma- 
tion of a short-lived and recurrent spiral structure, two episodes of 
which last over roughly 1 Gyr. In between these two episodes the 
spiral never quite vanishes, but its amplitude decreases consider- 
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ably. These spirals have the same properties as the manifold-spirals 
discussed in Papers I to V. We also followed the trajectories of par- 
ticles in the arms and found that, contrary to what would have been 
expected for density wave spirals, they do not cross the arms but 
they move along them, starting off from the vicinity of one of the 
Lagrangian saddle points {L\ or £2) in the direction of the other 
one. Thus their trajectories outline the arms. In fact the whole of 
the trajectory contributes to the spiral structure and these spirals 
can be called flux-tube manifold spirals. 

There are several more signs, tell-tale of a manifold origin. For 
example, particles which outline the spiral have their origin in the 
outer part of the bar and join the spiral via the vicinity of a point 
whose location is compatible with that of the saddle Lagrangian 
points (Li or L2). Furthermore, they loop in the vicinity of that 
point before they follow the arm. Another tell-tale sign of the man- 
ifold origin of the arm is that in the vicinity of each Lagrangian 
point the particles split into two groups, the same in simulations 
and in simple orbital calculations. 

From all the above it becomes clear that we are witnessing 
manifold-driven spirals in our simulation. This - together with the 
good agreement between manifold-driven and observed rings and 
spirals found in Papers IV and V and by Martmez-Garcta l |2012) 
- argues strongly that manifolds do play a role in spiral and ring 
formation. Manifolds, however, are not the only possible origin 
of such structures. Indeed, spirals have been witnessed also in 
other N-body simulations and interpreted in terms of other theo- 
ries (|Sellwood & Carlberg|1984"l|D'Onghia, Vogelsberger & Hern-| 
|quist|2012||Grand, Kawata & Cropper|2012||Serfwood|2012| >. We 
already included a note of caution to this avail in Paper V. 

Our comparison between N-body and theoretical manifolds 
must necessarily stay qualitative. As any orbital structure work, the 
manifold theory relies on a few simplifications, the most important 
of which is that the potential is time-independent in the frame of 
reference co-rotating with the bar. On the other hand, in simulations 
and in real galaxies the potential evolves with time, due to redistri- 
bution of angular momentum via the resonances (e.g. |Lynden-Bell| 
|& Kalnajs|1972)|Weinberg|1985HAthanassoula|2002|p003^ . How- 
ever, if this evolution is not too fast, it should not present a problem 
for our flux-tube manifold theory where the arm is constituted by 
the whole flow of material guided by the manifold^] and will only 
lead to a change of the manifold properties with time. Any secular 
change in the potential is expected to lead to a secular evolution of 
the manifold properties. Furthermore, our simulations show that, 
even when the rate of change is considerable, particle trajectories 
are still guided by manifolds, keep their characteristic signatures 
and obey the selection rules of permissible paths defined by them. 
Nevertheless, their shape and extent change as expected, while par- 
ticles can get trapped or untrapped by these manifolds. 

Such orbits as described here have also been witnessed as re- 
sponses to applied analytical potentials (e.g. [Danby 1965 , |Patsis| 
|2006| Papers III & V), although links to manifolds were not nec- 
essarily made. To our knowledge, however, this is the first time 
that the manifold theory is tested in a realistic self-consistent N- 
body simulation by following the motion of individual particle^] 



2 This is not the case for the theory presented in e.g. |Voglis, Tsoutsis 
& Efthymiopoulos ( 2008 1, or Harsoula, Kalapotharakos & Contopoulos 
(2011) which considers the loci of the apsidal manifolds sections. It thus 
requires a quasi-station ary, non-evolving potential. 

3 A previous attempt (Tsoutsis, Efthymiopoulos & Voglis|2008} used an 
unrealistic simulation, with no separate disc and halo components, just a 
single cylindrical-shaped component, with a vertical to horizontal size ratio 



Nevertheless, the simulation discussed here is in no way unique; 
manifold-driven spirals can be seen in a large number of our sim- 
ulations, some of which include gas, and will be discussed else- 
where. 

Comparing the disc extent between times 0.5 and 1.5 Gyrs we 
see a considerable change, of the order of 50%, due to the spirals. 
This shows that the bar and the associated manifolds can drive the 
overall evolution of the disc significantly. We will explore this pro- 
cess further in future work. 
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